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ABSTRACT 

Mechanical  admittance  methods  are  applied  to  two  simple 
composite  structures.  The  mechanical  admittances  of  the  com¬ 
posite  structures  are  predicted  using  the  mechanical  admittances 
measured  on  their  constituent  substructures.  The  effects  of  the 
suspension  system  on  the  measured  admittances  are  estimated. 
Attention  is  devoted  to  the  technique  of  analytically  assembling 
substructures  into  a  composite  structure. 

Resonant  frequencies  and  mode  shapes  of  the  composite 
structures  are  estimated  from  the  predicted  admittances.  The 
effects  of  modal  participation  and  damping  on  the  estimated  mode 
shapes  and  resonant  frequencies  are  discussed. 
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UHAPTER  I 
INTRODUCTION 


* 

The  mechanical  admittance  method  or  equivalently  the 
mechanical  impedance  method  is  a  means  of  cnaracterizing  the 
dynamic  behavior  of  a  structure  by  observing  its  dynamic 
behavior  at  a  finite  number  of  locations.  In  general,  admit¬ 
tance  and  impedance  are  complex  quanities  defined  for  steady 
state,  sinusoidal  motion  and  are  functions  of  the  frequency 
of  excitation.  The  admittance  method  Is  primarily  experi¬ 
mental  in  nature  and  versatile  in  its  applications. 

Mechanical  admittance  embraces  such  diverse  areas  as 
accoustical  transmission  in  structures  to  mode  shape  and 
resonant  frequency  determination.  Rubin'*’  uses  admittance 
methods  to  reduce  levels  of  engine  noise  transmitted  to  pas¬ 
senger  compartments  on  commercial  aircraft.  Admittance  meth¬ 
ods  are  also  applied  to  simple  structures  to  determine  reso- 

2 

nant  frequencies  and  transmission  matrices. 

The  fundamental  difference  between  an  admittance  mea¬ 
surement  and  an  impedance  measurement  is  the  displacement 
boundary  conditions  imposed  on  the  structure  during  the  mea¬ 
surement.  For  the  purposes  of  this  paper  impedances  are 
defined  as  the  elements  of  the  inverse  of  the  admittance 
matrix.  For  an  in  depth  treatment  of  the  differences  be¬ 
tween  admittance  and  impedance  see  reference  3- 


*  Admittance  is  synonymous  with  the  term  mobility. 
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This  paper  treats  the  problem  of  predicting  the  admit¬ 
tances  for  a  composite  structure  from  the  measured  admit¬ 
tances  of  the  constituent  substructures.  In  addition,  code 
shapes  and  resonant  frequencies  are  estimated  from  the  pre¬ 
dicted  admittances  of  the  composite  structure. 

Chanter  II  deals  with  the  mathematics  involved  in  the 
admittance  prediction  technique.  Chapter  III  is  a  brief  de¬ 
scription  of  the  total  experimental  system  and  the  function 
of  the  computer  programs.  Described  in  Chapter  IV  is  the 
use  of  the  total  system  for  predicting  the  admittances,  mode 
shapes  and  resonant  frequencies  for  two  simple  composite 
structures . 


ih 


CHAPTER  II 
ANALYTIC  BASIS 

2.1  The  Mechanical  Admittance  and  Mechanical  Impedance  of 

a  Structure 


For  most  structures  of  Interest  their  low  frequency 
dynamic  behavior  can  be  modeled  by  a  lumped  parameter  or 
a  finite  element  idealization.  Attention  is  focused  on 
linear  elastic  structures  undergoing  steady  state,  sinusoidal, 
forced  excitation.  This  approach  facilitates  the  derivation 
of  general  expressions  for  admittance  and  impedance. 


{qHm]{u}  +  [k]{u} 


Equation  2.1  is  the  general  equation  of  motion  for  a 
linear  elastic  structure.  (<5}  is  a  column  vector  of  the 
applied  forces  and  moments.  {u.}  is  a  column  vector  of  the 
resultant  local  displacements  and  rotations  of  the  structure. 
The  mass  matrix,  [M]  ,  and  the  stiffness  matrix,  [k]  ,  are 
square  symmetric  matrices  of  the  same  order  as  the  number  of 
degrees  of  freedom  of  the  structure. 


Schematic  of  the  Degrees  of  Freedom  of  the  Finite 


Element  Idealization  of  a  Continuous  Structure 

Figure  2.1 
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Equation  2.1  is  rewritten  for  steady  state,  sinusoidal 


motion. 


(o\  «  fa\eilrt 

{u}  -  {u}  ** 

(u1!  =  -  ul(u]  e1"* 


(2.2) 


(2.3) 


(2.4) 


(2.5) 


If  the  mass  and  stiffness  matrices  are  of  order  N,  then 
there  exist  N  nonogeneous  solutions,  (n  of  equation  2.2 


(o)  - 

.  {f}. 

ber  of  non-negative  frequencies  of  motion,  U)n.  Gow  is  de¬ 
fined  as  the  n—  resonant  frequency  of  the  structure.  me 


Tne  homogeneous  solutions 


(2.6) 


exist  only  for  a  finite  num- 


column  vector,  in  is  defined  as  the  n-^  mode  shape  of  the 
structure . 

In  general  the  mode  shapes  are  orthogonal  to  each  other 
with  respect  to  both  the  mass  and  stiffness  matrices.  Equa¬ 
tions  2.7  and  2.8  are  the  definitions  of  orthogonality. 


th 


LtfiWW  =  o. 

kwNM  -  O' 


J*K 


(2.7) 


(2.3) 


1 6 


The  resultant  motion  of  the  structure  for  forced  exci¬ 
tation  is  assumed  to  be  a  linear  combination  of  the  mode 
shapes . 

M  *  l2-9) 

The  elements  of  the  column  vector,  M.  are  the  modal  par¬ 
ticipation  functions  whicn  are  in  general  functions  of  fre¬ 
quency.  If  the  modal  participation  functions  are  determined, 
then  the  forced  response  of  the  structure  can  be  determined. 

Substituting  equation  2.9  into  equation  2.2  and  prenul- 
tiplying  both  sides  of  equation  2.2  by  W  ,  equation  2.10 
results . 

Invoking  tne  orthogonality  of  the  mode  shapes  with  respect 
to  the  mass  and  stiffness  matrices,  equations  2.7  and  2.8, 
the  mass  and  stiffness  matrices  are  transformed  into  diag¬ 
onal  matrices. 

[4>fM  =  +  [K]](a\  i2-11) 


Eacn  modal  participation  function  can  now  be  solved  for  al- 
ge  braically . 


+  k») 


(2.12) 


(2.13) 


(2.14) 
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Substituting  the  expression  i'or  the  modal  participation 
function,  Q.n  ,  into  equation  2.9  one  can  construct  the  forced 
response  of  the  structure. 


Um  -  ume 


LU>t 


=  4  L4>J  Wj  €elu>t 

^  (-oAn*  +  KyO 


(2.15) 


For  steady  atate,  sinusoidal  notion  the  resultant  velocity 

i.  ^ 

at  the  m~  degree  of  freedom  of  the  structure  ,"V_,  can  be 

9  w  It)  9 

expressed  as  a  function  of  the  resultant  displacement,  u  m » 
and  the  frequency  of  excitation,  U)  . 

Vw=  iu>umewt=  e  (m  k  I  2...n)  (2.16) 


(2.17) 


The  admittance  of  a  structure,  W  ,  is  defined  as  the 

ratio  of  the  resultant  velocity  at  the  m—  degree  of  freedom 

*1. 

of  the  structure  to  the  force  at  the  l-1-  degree  of  freedom, 
assuming  all  other  forces  equal  to  zero. 

_  L fO  Mix  6  \  TDy\  irn  T jL  .  ^  -  ,, » 

This  expression  for  admittance  is  valid  for  applied 
moments  and  rotational  velocities,  experimentally,  admit¬ 
tances  are  usually  measured  only  for  point  forces  and  trans- 


*  If  the  structure  is  modeled  as  a  continuum,  the  expres¬ 
sion  for  admittance  Is  altered  slightly.  is  replaced 

+ sf\  n  n 

by  the  n — -  continuous  mode  shape,  r  ,  evaluated  at  the  n — 

degree  of  freedom  of  the  structure. 
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lational  velocities.  For  the  remainder  of  the  paper  admit¬ 
tances  are  defined  for  point  forces  and  translational  veloc¬ 
ities  . 


The  admittances,  j[  (h>)  ,  are  usually  manipulated  in 
matrix  form,  since  admittances  are  measured  on  the  structure 
at  a  finite  number  of  points.  Using  superposition  of  veloc¬ 
ities,  an  admittance  matrix  relating  a  column  vector  of  point 
forces,  m  ,  to  a  column  vector  of  translational  velocities, 
,  can  be  formulated  as  in  equation  2.19. 


/  \ 
VI 

1  Vi 
< 

# 

^  ~ 

Y2l(u})  Ym(u>) 

.  TM 

<  •  > 
0 

• 

lL 

V 

YC“)  . 

XN| 

_ 

X 

(2.19) 


(2.20) 


All  of  the  admittances  must  be  evaluated  at  the  same 
frequency,  W  ,  for  the  relationship  to  hold.  Elements  on 
the  principal  diagonal  of  the  admittance  matrix  are  the  mu¬ 
tual  admittances.  All  other  elements  of  the  admittance  ma¬ 
trix  are  the  cross  admittances.  Using  equation  2.18  one  can 
show  that  the  admittance  matrix, 


Y 


.[Y]  ,  is  symmetric. 

N  ^  /bn  An  Jl-  cO  if*  in 

(.CO)  =  T  a  t  =  X  Tfrfrfrft  (2.21) 

The  p—  column  of  the  admittance  matrix, Y|ji*Hj  =  W)  , 
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can  be  interpreted  as  the  resultant  velocities  at  the  N  de- 
grees  of  freedom  for  a  unit  sinusoidal  force  at  the  p —  de¬ 
gree  of  freedom  with  no  forces  applied  at  the  other  degrees 

till 

of  freedom.  Deleting  the  p—  row  and  column  of  the  admit¬ 
tance  matrix  is  equivalent  to  setting  the  p—  force  equal  to 
zero . 


The  inverse  of  the  admittance  matrix  is  defined  as  the 
mechanical  impedance  matrix,  [2). 


M  -  [Ym  - 


Since  the  admittance  matrix  is  symmetric,  its  inverse,  the 
impedance  matrix,  is  also  symmetric.  The  elements  on  the 
principal  diagonal  of  the  impedance  matrix  are  the  mutual 
impedances.  All  other  elements  of  the  impedance  matrix  are 
the  cross  impedances. 


ft] 

4 

Z  „(<•>)  Zac«i  . 

.  .  Zm 

IN 

\ 

1 

{• 

4. 

« 

• 

Nl 

NH  -1 

i 

•  } 
t 

• 

irH 
{  * 

The  p~  column  of  the  impedance  matrix,  T—yp*  (j«  ^2....  bJ), 
can  be  interpreted  as  the  resultant  forces  at  the  N  degrees 

j.  y. 

of  freedom  for  a  unit  sinusoidal  displacement  at  the  p—  de¬ 
gree  of  freedom  with  no  velocity  allowed  at  the  other  degrees 
of  freedom.  Deleting  the  p—  row  and  column  of  the  impedance 
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matrix  Is  equivalent  to  setting  the  p—  velocity  equal  to 
zero . 

Prom  a  practical  point  of  view,  making  admittance  mea¬ 
surements  on  a  structure  is  easier  than  making  impedance 
measurements.  To  make  an  admittance  measurement  a  force  is 
applied  at  one  degree  of  freedom  and  the  resultant  velocity 
is  monitored  at  a  degree  of  freedom  with  no  velocity  (dis¬ 
placement)  boundary  conditions  imposed  on  the  structure.  To 
make  an  impedance  measurement  a  velocity  must  be  produced  at 
one  degree  of  freedom  while  the  remaining  degrees  of  freedom 
are  held  fixed.  The  reaction  forces  at  the  fixed  degrees  of 
freedom  as  well  as  the  force  at  the  point  of  excitation  must 
be  measured. 

The  stiffness  matrix,  M  ,  and  the  flexibility  matrix, 
[cl  ,  used  for  static  analysis  of  structures,  are  analogous 
to  the  impedance  matrix,  tel  ,  and  the  admittance  matrix, 

[Y]  ,  respectively.  Compare  equations  2.22  and  2.20  with 
equations  2.24  and  2.25. 


([V  MM 

M  is  a  column  vector  of  static  loads  and 
vector  of  static  displacements. 


(**}  i 


(2.24) 


(2.25) 


s  a  column 
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2,2  Assembling  Substructures  Into  a  Composite  Structure 

Using  Mechanical  Impedance  Notation 

Assuming  that  admittance  matrices  have  been  measured 
for  a  number  of  substructures  to  be  physically  assembled 
Into  a  composite  structute,  one  can  predict  the  properties 
of  the  composite  structure  by  mathematically  assembling  the 
substructures  using  Impedance  notation.  If  the  substructure 
admittance  matrices  are  Inverted  to  substructure  Impedance 
matrices,  the  substructure  Impedance  matrices  can  be  added 
to  construct  the  composite  structure  Impedance  matrix. 


Region  d^  Region  c ^  Region  02  Region  d2 


Substructure  #1  Substructure  #2 

Schematic  of  Substructures  to  be 
Assembled  Into  a  Composite  Structure 

Figure  2.2 

Each  substructure  In  Figure  2.2  Is  partitioned  by  a 
broken  line.  Those  degrees  of  freedom  that  are  to  be  Joined 
to  degrees  of  freedom  on  the  other  substructure  are  parti¬ 
tioned  from  the  other  degrees  of  freedom  in  the  substructure. 

The  substructure  impedance  matrices  in  equations  2.26 
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and  2.27  are  partitioned  in  an  analogous  manner.  The 
elements  of  each  impedance  matrix  corresponding  to  degrees 
of  freedom  that  are  to  be  Joined  to  degrees  of  freedom  on 
the  other  substructure  are  partitioned  from  the  elements  of 
the  matrix  corresponding  to  the  other  degrees  of  freedom  of 
the  substructure.  Regions  c^  and  C2  In  Figure  2.2  must  have 
the  same  number  of  degrees  of  freedom  to  insure  that  f 

all  have  the  same  number  of  elements. 


KtV  and  {4^ 
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.  2C, d,  ;2Vi 
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(2.26) 


f, 


C2I 


4 


cz°i 


doC 


r-z 


Z 


cz32 


d2d 


203J 


Vfc, 


(2. 27) 


Equations  2.26  and  2.27  are  separated  into  four  matrix 


equations  for  the  forces  in  the 

four  regions 

of  Figure  2.2. 

- 

(2.28) 

II 

oTc,} 

(2.29) 

(U  • 

( z<,  - 

(2.30) 

{•Q  = 

( z,,,,]  w 

+ 

)W 

(2.31) 
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Using  equilibrium  between  Internal  forces,  (Q.m 
and  applied  forces,  (q.  equations  2.29  and  2.30  can  be 
added  to  give  the  applied  forces  In  region  c  of  the  composite 
structure. 


{Q  -  (Q  =  {*} 


(2.32) 


(2.33) 

+  [z  czdzlfe} 


Region  d^ 


Region  c 


Region  d2 


Schematic  of  Assembled  Composite  Structure 

Figure  2.3 

Invoking  compatibility  of  velocities,  equation  2.34,  at 
the  degrees  of  freedom  common  to  both  substructures  in  region 
c,  equation  2.33  can  be  reduced  to  equation  2.35. 

{ifc,}  =  { * {V}  (2-3+) 

M  =  +  [Zc,c,+  Zc2c2V^  +[  ZcjdilM  <2.35) 


The  impedance  matrix  for  the  composite  structure  is 


2k 


constructed  from  equations  2.28,  2.31  and  2.35. 
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(2.36) 


Since  the  measured  admittances  for  the  substructures 
are  valid  only  for  translational  velocities,  the  substructure 
impedance  matrices  are  valid  only  for  translational  velocities. 
Consequently,  at  the  points  of  attachment  of  the  substruc¬ 
tures  one  can  insure  the  compatibility  of  translational 
velocities  only,  equation  2.34.  It  follows  that  the  attach¬ 
ment  points  of  the  first  substructure  are  allowed  to  rotate 
Independently  from  the  attachment  points  of  the  second  sub¬ 
structure.  As  a  result  of  this  limitation  in  assembling 
substructures  analytically,  composite  structures  that  have 
large  rotational  impedances  at  the  attachment  points  of 
their  constituent  subsrtuctures  will  not  be  accurately 
modeled  by  the  composite  structure  impedance  matrix. 

Displacement  boundary  conditions  are  imposed  on  the 

composite  structure  by  operations  on  the  composite  struc- 

ture's  Impedance  matrix.  One  makes  the  velocity  of  the  p — 

til 

degree  of  freedom  zero  by  deleting  the  p —  row  and  column 
of  the  Impedance  matrix.  Since  impedances  are  available 
Just  for  translational  velocities,  only  translational 
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velocities  of  the  composite  structure  can  he  made  zero. 

After  the  displacement  boundary  conditions  are  Imposed 
on  the  composite  structure  Impedance  matrix,  it  is  Inverted 
to  the  composite  structure  admittance  matrix.  The  resonant 
frequencies  and  accompanying  mode  shapes  of  the  composite 
structure  are  estimated  from  the  predicted  admittance  matrix 
for  the  composite  structure  as  shown  in  section  2.4. 
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2,3  The  Effect  of  Damping  on  Resonant  Frequencies  and 

Admittances 


In  section  2,1  the  resonant  frequencies  of  a  structure 
are  defined  as  tOn  and  the  mode  shapes  are  defined  as 
As  the  frequency  of  excitation,  CO  ,  tends  to  GOn  ,  the 
forced  response  of  the  structure  tends  to  Infinity  where  the 
ratio  of  the  response  at  one  degree  of  freedom  of  the  struc¬ 
ture  to  the  response  at  another  degree  of  freedom  Is  defined 
by  {4>%  The  property  of  damping  exhibited  In  most  real 
structures  changes  the  definition  of  resonant  frequencies. 
Consequently,  damping  affects  the  magnitude  and  phase  of  an 
admittance  at  or  near  a  resonant  frequency  of  the  structure. 

To  facilitate  the  discussion  of  damping  It  Is  assumed 
that  the  damping  present  In  the  structure  Is  modal  damping. 
Each  mode  Is  assumed  to  exhibit  damping  Independently  of  the 
damping  of  other  modes. 

This  form  of  damping  Is  represented  mathematically  by 
Introducing  a  diagonal  damping  matrix  Into  equation  2.11. 

hi  eiwt  =  [  -  [Ml  +  iu>te>3  +  [K]](al  e,ut  <a.37) 

The  modal  participation  functions,  Cln  ,  are  solved  for  in  the 
same  manner  as  in  section  2.1. 


(2.38) 
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«  ,  V  .  SStM 

^nl(o)  -  rrinTi 


(2.39) 


is  the  modal  damping  coefficient.  For  viscous 
damping  is  d„u>n  .  For  hysteretic  damping  is  of 

the  form  shown  in  equation  2.40. 

2 


*>>  =  ^ 


(2.40) 


The  admittance  of  a  structure  with  modal  damping  is 
given  in  equation  2.41. 


N  U>  A"  a.* 

YM=ym*M _ =y: 


h* 


-ull 


2.41) 


The  admittance  matrix  for  a  structure  with  modal  damping  is 
symmetric.  Consequently,  the  impedance  matrix  for  a  struc¬ 
ture  with  modal  damping  is  also  symmetric. 

The  poles  of  equation  2.41  are  the  damped  resonant 
frequencies  of  the  structure ,  51  n  ,  and  are,  in  general, 
complex  numbers. 


‘O# „)Sl„  +  t  (0*-  <*>* )  -  o. 


The  damped  resonant  frequencies  can  be  solved  for  explicitly 
for  viscous  damping,  equation  2.43,  and  hysteretic  damping, 
equation  2.44. 


-  to* 

(.  ■  d- 

U  "T/  +  “o*  V 

=  d„coK 

(2.43) 

Q-y\  ~  ( 

■  ■  .  x  H 

1  +  idn) 

(2.44) 
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The  damping  coefficient,  dn  »  is  much  less  than  one  for  low 
frequencies  of  excitation  in  most  structures.  For  hysteretic 
and  viscous  damping  is  the  q;  factor  or  the  resonant 

a.  u 

quality  of  the  n—  mode.  For  the  remainder  of  the  chapter 
q;  is  assumed  to  be  50. 

j.  y. 

If  the  frequency  of  excitation  is  near  the  J —  undamped 
resonant  frequency,  Wj  #  the  J —  term  in  equation  2.41  is 
most  representative  of  the  structure's  admittance. 


X'?>  ■ 


 m]  Jits. 


+  e 


tnSL 


(2.45) 


^ —  jyL  A*  U. 

r  a  9<n  ££ _ 

tnJL  nTT  + 


(2.46) 


£mjL  represents  the  participation  of  all  modes  but  the 
J —  for  .  The  admittance  matrix  for  the  structure 

excited  near  the  J —  undamped  resonant  frequency  is  given  in 
equation  2.47. 
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(2.48) 


If  all  are  zero  for  CO»tOj  ,  then  the  matrix,  IT), 

is  equal  to  the  admittance  matrix.  All  of  the  elements  of 
(t]  are  either  In  phase,  0  radians,  or  out  of  phase,  TC 
radians,  with  respect  to  each  other  for  CO » (*)'  .  All  the 
elements  of  fT)  have  a  local  maximum  at  one  frequency  of 
excitation.  This  frequency  Is  the  real  part  of  the  damped 
resonant  frequency  defined  In  equation  2.42. 

In  general  the  participation  of  all  modes  except  the 
J— ,  is  n°t  zero  foru)»U)j.  It  follows  that  the  ele¬ 
ments  of  the  admittance  matrix,  +  *  w111  not  have 

their  local  maxima  at  exactly  the  same  frequency  for  . 

In  addition  the  elements  of  the  admittance  matrix  are  not 
exactly  in  or  out  of  phase  with  respect  to  each  other  for 
any  frequency  of  excitation. 
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2,4  Estimating  Resonant  Frequencies  and  Mode  Shapes  for 

Composite  Structures  from  Predicted  Admittance  Data 


The  predicted  admittances  of  a  structure  are  the  normal¬ 
ized  levels  of  the  forced  response.  The  levels  of  the  forced 
response  are  normalized  by  the  magnitudes  of  the  applied 
forces. 

Equation  2.47  Is  representative  of  the  behavior  of  an 

th 

admittance  matrix  at  or  near  the  J —  undamped  resonant  fre¬ 
quency  of  the  structure. 
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(2.48) 


til 

If  the  participation  of  all  modes  except  the  J — ,  , 

X  L- 

Is  very  small  compared  to  the  participation  of  the  J —  mode, 

x.  y. 

then  one  can  predict  the  resonant  frequency  and  mode 
shape  from  the  admittance  matrix.  The  frequency  at  which 
the  elements  of  the  admittance  matrix  have  their  J —  local 
maximum  Is  the  resonant  frequency.  The  rows  and  columns 
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of  the  admittance  matrix  are  proportional  to  the  J —  mode 
shape,  W).  for  ft:  u>j  . 

In  general  there  are  two  cases  where  the  participation 

til 

of  all  modes  except  the  Is  not  small  compared  to  the 

Vi 

participation  of  the  J —  mode  for  (Ocs  u>j  .  First,  if  an 
admittance  Is  predicted  for  a  degree  of  freedom  at  or  near 

j.  u 

a  node  of  the  mode  shape,  that  admittance  Is  more 

representative  of  the  participation  of  modes  other  than  the 
th 

3 — .  Second,  if  two  resonant  frequencies  of  the  structure 
are  close  together,  the  admittance  matrix  is  representative 
of  the  participation  of  both  modes  at  either  resonant 
frequency. 

If  an  element  of  the  admittance  matrix  corresponds  to  a 
degree  of  freedom  near  a  node  of  the  mode  snape  and  Is 

i.  y. 

dominated  by  the  participation  of  modes  other  than  the  J — , 
then  that  element  will  not  nave  a  well  defined  maximum  in 


magnitude  in  the  region  of  the  J —  resonant  frequency. 


N 

The  expression  for  the  participation  of  all  modes  but  the 

"t  ^  O 

cm&  ,  does  not  have  any  poles  near  U>j  .  Therefore , there 
is  no  local  maximum  In  magnitude  of  £mj.  for  . 

If  two  resonant  frequencies  of  the  structure  are  close 
together,  then  the  admittance  can  be  approximated  by  the 
participation  of  only  tnose  two  modes  for  frequencies  of 
excitation  near  the  two  resonant  frequencies. 
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(2.49) 


Evaluating  the  admittance  at  GOj  one  can  solve  for 
U3p 

values  of  ttt  such  that  the  magnitude  of  the  participation 

, ,  J  +h 

of  the  p—  mode  is  oC  of  the  participation  of  the  J-=—  mode 

th  „ 

at  the  J —  resonant  frequency. 

oC 


$4>p  <t>[ 

Wlp  Tni  rJL 

%(u>^  i(co?-  u> p1) 

(2.50) 


The  damping  for  both  modes  is  assumed  to  be  viscous. 

%(ui)  =  djIOj  =  Zr^Oj  (2.5D 

=  dplOp  =  2y^,Up  (2.52) 

is  the  ratio  of  the  damping  of  the  mode  to  critical 
damping  for  that  mode.  The  ratio  of  the  normalized  modal 
products  for  the  p—  and  modes  is  |3  • 
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-STj  (ft, 


(2.53) 


Substituting  equations  2.51»  2.52  and  2.53  into  equa- 

CJj 
U>\ 

1  rX 

(2.54) 


tion  2.50  one  can  solve  for  values  of  — c 


X*-  (2-^")X  +  1  -  4^  =0 
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A  =  oo.2- 
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(2.55) 


33 


Plots  of  Equation  2.54-  for  Different 
Values  of  the  Parameters  dn  ,  oc  »  andjS 

Figure  2.4 
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Figure  2.4  shows  plots  of  equation  2.54  for  different 
values  of  the  parameters,  »  O C  and  .  Assuming  the 
modal  damping  coefficients  dn  are  ,02,  the  ratio  of  the 
normalized  modal  products  ^  is  1.0  and  0(.  is  .1,  the 
roots  of  equation  2.54  are  given  in  equations  2.56  and  2.57. 

cof 

X  =  <*>}*  ~  1.2  (2.56) 
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(2.57) 


Thus, if  is  greater  than  .8  and  less  than  1.2,  the 
participation  of  the  p~  mode  at  the  J—  resonant  frequency 

tOf?" 

is  greater  than  10%.  If  TTi  is  less  than  .8  or  greater  than 

j 

xu  xu 

1.2,  the  participation  of  the  p—  mode  at  the  J—  resonant 
frequency  is  less  than  10%. 

The  participation  of  the  p~  mode  affects  the  phase  of 

^u 

the  admittance  at  the  J —  resonant  frequency.  The  two  modal 
participation  terms  in  equation  2.49  are  treated  as  vectors 
in  the  complex  plane. 

%  € 
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(2.49) 


The  Modal  Participation  Vectors  of  the  P-^=  and  J~  Modes  in 


the  Complex  Plane 
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The  phase®  of  the  two  modal  participation  vectors  evaluated 
at  (Oj  are  given  in  equations  2.58  through  2.61. 


=  -ton''  =0.  (jt  <14  >  O 

=  TT  -  tan'  =  TC  M< 0 

1  $A>0 


(2.58) 

(2.59) 

(2.60) 

(2.61) 

If  is  the  magnitude  of  the  participation  of  the 

mode  and  M  is  the  magnitude  of  the  participation  of  the  p— 

^  th 

mode,  the  phase  of  the  admittance  at  the  J— -  resonant  fre- 

quency  is  given  in  equations  2.62  and  2.63. 
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(2.62) 


(2.63) 


The  participation  of  the  p~  mode  affects  the  phase  of 


the  admittance  most  when  ([)  Vi>=±¥  . 


(2.64) 


(j)  =  tan"'  -  tan”'|i(x|  <j)J(^j)  =  0 

(E>(a>i)  =  TC  +tan"'  j±j^l  =  tan  V+otv  &TC±<x  .  (2.65) 

1  niJ  1 


36 


th 

The  effect  of  the  p —  mode  on  the  phase  of  the  admittance  at 

the  ^ —  resonant  frequency  is  at  most  OC  radians.  OC  Is 

defined  In  equation  2.50. 

th 

Assuming  that  the  J—  resonant  frequency  Is  not  close 

Vi 

to  other  resonant  frequencies  of  the  structure,  the  J —  mode 
shape  Is  calculated  from  the  maximum  values  of  the  magnitudes 

‘V  Vi 

for  the  admittances  near  the  3 —  resonant  frequency.  An  ele¬ 
ment  of  the  admittance  matrix  that  does  not  have  a  local 
maximum  In  magnitude  for  and  Is  dominated  by  the 

participation  of  modes  other  than  the  J —  is  not  used  to 

th 

calculate  the  J—“  mode  shape. 

The  signs  of  the  components  of  the  mode  shape  vector, 

are  determined  by  the  relative  phases  of  the  admit- 
th 

tances  at  the  J —  resonant  frequency.  Unfortunately,  the 

admittance  data  is  available  only  for  a  finite  number  of 

frequencies  which  may  not  coincide  with  the  exact  resonant 

frequencies  of  the  composite  structure.  However,  as  long  as 

th 

the  participation  of  modes  other  than  the  J —  is  small  for 

th 

LOSsWj  ,  the  relative  phases  of  the  admittances  at  the  j — 

resonant  frequency  are  approximately  the  same  as  their 

th 

relative  phases  near  the  J —  resonant  frequency. 


The  matrix 


[J] 


is  approximately  equal  to  the  admlt- 


.th 


tance  matrix  at  the  J —  resonant  frequency  if  the  partlcl- 

th  # 

patlon  of  modes  other  than  the  J —  Is  small.  Equations  2 .( 
and  2.67  give  the  phases  of  elements  of  the  matri::  CU  . 


*  See  equation  2.48 
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€  £>  o 

(2.66) 

££<o 

(2.67) 

are 


The  relative  phases  of  the  elements  of  the  matrix  m 
0  or  TC  radians  for  all  frequencies  of  excitation, CO  .  These 
relative  phases  determine  the  signs  of  the  components  of  the 
mode  shape  vector,  and  are  approximately  equal  to 


the  relative  phases  of  the  admittances  at  the  resonant 
frequency  of  the  structure. 
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CHAPTER  III 


The  Mechanical  Admittance  Measurement  and  Prediction  System 

The  primary  function  of  the  measurement  and  prediction 
system  is  to  predict  the  admittances  of  a  composite  structure 
from  admittance  measurements  made  on  the  constituent  sub¬ 
structures.  The  admittance  measurements  for  the  substruc¬ 
tures  are  made  in  analog  form  and  then  converted  to  digital 
form  to  construct  the  substructure  admittance  matrices. 

The  substructure  admittance  matrices  are  used  to  calculate 
the  composite  structure  admittance  matrix. 

The  admittances  of  a  substructure  are  constructed  from 
force  and  acceleration  measurements  made  on  the  substruc¬ 
ture.  The  force  and  acceleration  measured  are  resultants 
of  a  force  input  to  the  structure  generated  by  an  electro¬ 
magnetic  shaker.  The  resultant  force  imparted  to  the  sub¬ 
structure  is  measured  by  a  very  stiff,  lightweight  strain 
gauge  mounted  between  the  shaker  and  the  substructure.  Tne 
resultant  acceleration  is  measured  by  an  accelerometer 
mounted  behind  the  strain  gauge  or  at  some  location  on  the 
substructure. 

The  force  input  to  the  structure  generated  by  the  shaker 
is  proportional  to  the  deflection  of  the  strain  gauge. 

Because  the  alloy  in  which  the  strain  gauge  is  embedded  is 
stiff  and  lightweight,  the  acceleration  measured  benind  the 
strain  gauge  is  the  local  acceleration  of  the  substructure 
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The  Mechanical  Admittance  Measurement  and  Prediction  System 

Figure  3.1 


4o 


in  front  of  the  strain  gauge  ezcept  at  very  high  frequencies 
4 

of  excitation. 


ACCELERATION  Output  (A) 


Schematic  of  Shaker  and  Impedance  Head 
Figure  3.2 

A  strain  gauge  with  an  accelerometer  mounted  behind  it 
as  in  Figure  3.2  is  referred  to  as  an  "impedance  head".  The 
device  used  in  the  following  experiments  is  a  Z602  impedance 

head  made  by  Wilcoxon  Research.  The  shaker  is  a  F3  "wrap 

* 

around”  shaker  also  made  by  Wilcoxon  Research. 

The  admittances  of  a  substructure  are  measured  over  a 
finite  frequency  range  by  varying  the  frequency  of  the  input 
to  the  shaker  from  the  SD104A  sweep  oscillator.  The  resul¬ 
tant  force  and  acceleration  signals  are  amplified  and  put 
into  the  analog  impedance  equipment.  These  signals  are  band 
pass  filtered  before  they  are  used  to  calculate  the 

log  (magnitude)  and  phase  of  the  measured  admittance. 

10 

Each  SD101A  filter  is  made  of  many  band  pass  filters. 


*  The  specifications  for  these  instruments  appear  in 
reference  5. 


41 


Schematic  of  the  SD1002B-33  Automatic  Mechanical  Impedance 
System  In  the  Mechanical  Admittance  Mode 

Figure  3.3 


OUTPUTS  To 
S-WAt-Olj  r'LO.TEK 


*  For  a  complete  description  of  the  SD1002B-33  Automatic 
Mechanical  Impedance  System,  see  reference  7. 
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Each  band  pass  filter  has  a  unique  center  frequency  and  is 
active  over  a  finite  frequency  range  that  Includes  its 
center  frequency.  The  active  ranges  of  the  band  pass  fil¬ 
ters  cover  a  broad  frequency  range,  5  Hz.  to  10  kllo-Hz., 
and  they  do  not  overlap. 


ACTIVE  RavG£  Active  Kange 
of  fiRST  filter  of  second  filter 


i.  - 


Gain 


iotc 

SEC. 


FREQUENCE  of 

Excitation 


Schematic  of  the  Transfer  Function  of  a 
SD101A  Filter  in  the  Frequency  Domain 
Figure  3.4 

Only  one  band  pass  filter  is  active  for  a  given  fre¬ 
quency  of  excitation  generated  by  the  SD104A  sweep  oscil¬ 
lator.  The  resultant  force  and  acceleration  signals 
generated  by  the  impedance  head  should  have  the  same  fre¬ 
quency  as  the  frequency  of  excitation  for  the  shaker 
generated  by  the  SD104A.  If  other  frequencies  are  present 
in  the  force  and  acceleration  signals  outside  the  band  of 
the  active  filter,  they  are  eliminated  from  the  signals. 
This  filtering  is  done  to  make  more  accurate  phase  and  mag¬ 
nitude  calculations  for  the  admittance. 


The  filtered  acceleration  signal  is  divided  by  the 
frequency  of  excitation  and  shifted  in  phase  by  -y  radians 
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to  yield  the  velocity  of  the  substructure.  The  velocity 
signal  and  the  filtered  force  signal  are  used  to  calculate 
the  admittance  of  the  structure.  The  SD112  log  converter 
generates  a  voltage  proportional  to  the  log^( magnitude )  of 
the  admittance.  The  329BSD  phase  meter  generates  a  voltage 
proportional  to  the  phase  of  the  admittance. 

The  voltages  proportional  to  phase  and  log^^magnitude) 
are  plotted  on  a  139AR  Hewlett  Packard  X  Y^  Y^  Recorder. 

These  voltages  are  also  preamplified  and  put  into  a  Radiation 
analog- to-digital  (A/D)  converter. 

The  A/D  converter  has  one  analog  input  channel  and 
generates  12  bit  digital  words  from  the  input  voltage.  Tne 
input  voltage  range  is  -2  volts  to +2  volts.  With  a  twelve 
bit  register  the  A/D  converter  can  generate  4099  different 
numbers  corresponding  to  the  four  volt  input  range  or  roughly 
one  number  for  every  millivolt  input. 

The  outputs  from  the  analog  impedance  system  are 
typically  50  millivolts.  To  achieve  any  accuracy  In  the 
digital  data  these  signals  must  be  preamplified  before  they 
are  put  into  the  A/D  converter.  After  making  a  preliminary 
frequency  sweep  to  establish  the  maximum  and  minimum  values 
of  the  voltage  output  from  the  Impedance  system,  the  bias 
and  gain  are  adjusted  on  the  preamplifier  to  take  full 
advantage  of  the  A/D  converter's  four  volt  input  range. 


*  For  a  complete  description  of  the  Radiation  A/D  converter 
see  reference  6. 
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The  SD1002B-33  Automatic  Mechanical 
Impedance  System  and  Charge  Amplifiers 


Figure  3.5 
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A  second  frequency  sweep  is  made  to  sample  and  convert  the 
amplified  analog  signal  from  the  impedance  system,  and  the 
digital  data  is  stored  on  magnetic  tape. 

Because  the  A/D  converter  has  only  one  input  channel 

the  log^Cmagnitude )  and  the  phase  of  the  admittance  are 

# 

converted  on  two  seperate  frequency  sweeps.  The  digital 
data  is  converted  in  format  from  12  bit  words  to  32  bit 
words  that  the  IBM  360  computer  can  tinders tand.  Each  12 
bit  word  is  padded  with  20  extra  bits  without  changing  the 
value  of  the  number  in  the  word  by  the  computer  program 
CONVERT. 

The  digital  data  is  scaled  from  the  analog  plots  made 
on  the  Hewlett  Packard  recorder.  Two  values  of  the  pnase  or 
magnitude  from  a  plot,  P-j^  and  P2,  are  compared  to  two  values 
of  the  digital  data,  and  D 2,  for  corresponding  frequencies 
of  excitation.  The  gain  and  bias  introduced  into  the  digital 
data  by  the  preamplifier  are  calculated. 


Gain  =  (P2-  P,) 


P,(Pi-P,) 

Blfts=  (Di-D,) 


(5.1) 

(3.2) 


*  If  a  two  track  tape  recorder  is  available  it  is  better  to 
record  phase  and  log^Q(magnitude )  on  the  same  frequency 
sweep  and  then  to  put  the  analog  signals  into  the  A/D  con¬ 
verter  one  track  at  a  time. 


b6 


All  of  the  digital  data  corresponding  to  the  analog  plot  is 
scaled  using  the  calculated  gain  and  bias.  D  is  the 
unsealed  value  of  the  digital  data  for  the  frequency  of 
excitation.  Sm  is  the  scaled  value  of  Dm. 

SM=  +  (Gmn)~'Dm  =  F»  -  (3.3) 

This  scaling  process  must  he  done  twice  for  each  admittance, 
once  for  phase  and  once  for  magnitude.  All  scaling  is  done 
by  the  computer  program  SCALE. 

The  scaled  data  is  used  to  construct  the  substructure 
admittance  matrix  for  the  sampled  frequencies  of  excitation. 

After  all  of  the  substructure  admittance  matrices  are  con- 

* 

structed  they  are  inverted  to  substructure  impedance  matri¬ 
ces.  The  substructure  impedance  matrices  are  assembled  into 
the  composite  structure  impedance  matrix.  All  matrix  inver¬ 
sion  and  assembling  is  done  by  the  computer  program  ADMIT. 

The  mode  shapes  and  resonant  frequencies  of  the  com¬ 
posite  structure  are  estimated  from  the  admittances  of  the 
composite  structure  as  shown  in  section  2.4. 


*  The  method  of  matrix  inversion  is  described  in  Appendix  A 
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CHAPTER  IV 
EXPERIMENTAL  RESULTS 

4.1  Introduction  to  Experimental  Results 


The  admittance  measurement  and  prediction  system  is 
applied  to  two  composite  structures.  The  first  composite 
structure  is  an  80-inch  free  rod  assembled  from  two  40-inch 
free  rods.  The  second  composite  structure  is  two  32-inch 
beams  pinned  together. 

32-inck  Beam 

© 

42-inch  Composite 
8eam 


32- inch  66AM 


Schematic  of  Assembling  Two  Composite  Structures 
from  Their  Constituent  Substructures 
Figure  4.1 

Admittance  measurements  are  made  only  on  the  substruc¬ 
tures.  The  "predicted"  admittances  are  the  admittances  of 
the  composite  structure  predicted  by  the  computer  program 
ADMIT.  The  "calculated"  admittances  are  generated  by 
analytical  modal  models  of  the  structures  and  are  used  only 
for  comparison  with  the  "measured"  and  "predicted"  admittances. 


Measuring  the  Cross  Admittance 
of  a  40-Inch  Steel  Rod 


Figure  4.2 
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4.2  The  Effects  of  the  Suspension  System,  Impedance  Head 

and  Shaker  on  the  Measured  Admittances  of  a  40-Inch  Free  Rod 


A  rubber  band  suspension  system  shown  in  Figure  4.2  is 
used  to  simulate  a  free  rod  for  frequencies  of  excitation, cO  , 
well  above  the  resonant  frequency  of  the  suspension  system. 
The  resonant  frequency  of  the  rod  and  rubber  band  ensemble 
is  1.2  c.p.s.,  which  is  well  below  the  lowest  frequency  for 
which  admittance  measurements  are  made  on  the  rod 

(2TC^  50 

c.p.s. ) . 


In  order  to  estimate  the  effect  of  the  rubber  band  on 
the  measured  admittances,  the  rod  and  rubber  band  are  modeled 
as  a  two  degree  of  freedom,  lumped  parameter  system.  See 
Figure  4.3.  Because  the  model  has  only  two  degrees  of  free¬ 
dom  it  will  be  representative  of  the  system  for  frequencies 
of  excitation  up  to,  and  including,  the  second  resonant 
frequency  of  the  rod  and  rubber  band  ensemble.  However, 
this  two-degree  of  freedom  model  will  give  some  Indication 
of  the  effects  of  the  rubber  band  suspension  system  on  the 
measured  admittances  of  the  rod. 


R 

itf 


"8 


kg  -  the  stiffness  of  the  rubber  band 


Y Y\ 


CB  “ 


m 


'  damping  " 
k^  -  the  stiffness  of  the  rod 
m  -  the  lumped  masses  of  the  rod 


The  Lumped  Parameter  Model  of  the  Rubber 
Band  and  Rod  Ensemble 


Figure  4.3 
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The  impedance  matrix,  [2]  ,  of  the  ensemble  is  written 
as  the  sum  of  a  mass  matrix,  [Ml  ,  a  damping  matrix,  !c}  , 
and  a  stiffness  matrix,  tK]  .  The  damping  in  the  rod  is 
assumed  to  be  zero;  the  damping  in  the  rubber  band  is  assumed 
to  be  viscous. 


[Z]  =  +[c)  +(i*of[K.] 


w  0. 

+ 

o.  o. 

Kr  -Kr 

_  o.  m  _ 

o.  c8_ 

;Kr  Kr*Xs 

The  admittance  matrix  for  the  system  is  found  by  invert¬ 
ing  the  impedance  matrix.  The  elements  of  the  admittance 
matrix  are  shown  in  equations  4.5  through  4.7,  and  the 
parameters  and  w2  are  defined  in  equations  4.3  and  4.4. 


W 


W2  = 


(4.3) 

(4.4) 


Y.w- 
Y„M  - 

Y.<">  - 
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)dL  -  (4.8) 

W2Z 

For  comparison  the  admittances  of  a  free  rod,  modeled 
by  a  spring  and  two  masses,  are  given  in  equations  4.9  and 
4.10. 

4 

<R 

f. 


Lumped  Parameter  Model  of  a  Free  Rod 
Figure  4.4 


<o  (u£  _  j_\ 

Y((o)  =  y  (to)  =  w  2  L 


(4.9) 


(4.10) 


is  the  first  resonant  frequency  of  the  free  rod 
for  elastic  motion. 


GO*  =  W2 


■m 


& 


(4.11) 


w,1 


If  and  Cg  are  zero,  the  admittances  for  the  rod 

and  rubber  band  ensemble  are  equal  to  the  admittances  for 

W,z 

the  free  rod.  For  the  system  shown  in  Figure  4.2,  r^a.  is 

5.6  X  10‘7  and  cp  is  .405  ■  lb;~3ec>  .  This  value  of  c^  is 

a  in.  B 

17$  of  critical  damping. 
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wt 

•^>2  ls  derived  using  a  calculated  stiffness  for  the  rod 
and  a  measured  stiffness  for  the  rubber  band. 


\A1  X  I0r 

lYV 


(4.12) 


kB  =  .33 


lb. 

)n. 


(4.13) 


W?  _  _Kb_ 

w2*  '  4Kr 


|0"7 


(4.14) 


cB  ls  calculated  by  observing  the  rate  of  decay  of  motion 
for  the  rod  and  rubber  band  ensemble  oscillating  In  its 
first  mode. 


W,* 

By  substituting  the  values  of  and  c^  Into  equations 
4.5  through  4.7  the  effect  of  the  rubber  band  support  system 
on  the  admittances  is  determined. 


Ym  = 

-Hi 


£i[w> -Kt *  405 


(4.15) 


=Yj“>  =  l(-  ,o-7)-W,‘)  +  (^Y¥Yo} 


(4.16) 


Y^>  = 


_  (4-17) 

'{%- <**-  >0'7)- W,2)  +  (^-  £)&  .403 


If  C*3  ,  the  circular  frequency  of  excitation,  is  much 


greater  than  w^,  the  admittances  are  affected  very  slightly 
by  the  stiffness  of  the  rubber  band.  At  the  first  resonant 
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frequency  of  the  free  rod,  the  magnitudes  of  the  admit¬ 

tances  are  very  sensitive  to  the  values  of  Cg. 

The  additional  mass  of  the  shaker  and  impedance  head 
attached  to  the  lower  end  of  the  rod,  Figure  4.2,  also 
affects  the  measured  admittances.  The  impedance  matrix  for 
a  two-degree  of  freedom, free  rod  with  the  mass  of  the  shaker 
and  impedance  head,  ms,  added  is  given  in  equation  4.18. 


[Z]=(i«) 


m  +  ms  q 
O.  m 


+  (to)-’ 


Kr 

~Kr 


(4.18) 


The  expression  for  the  first  resonant  frequency  for 


elastic  motion  is  CO 


=  00, 


Vz 


(4.19) 


If  1,  the  resonant  frequency  of  the  rod  with  the  shaker 

and  impedance  head  attached  is  not  appreciably  affected  by 
the  additional  mass,  mB.  However  ,  — *  for  the  ensemble  in 
Figure  4.2  is  .786,  and  the  first  resonant  frequency  for  the 
40-inch  rod  with  the  shaker  and  impedance  head  attached 
should  be  significantly  lower  than  CO* . 


00,  =  .88600* 


(4.20) 


By  modeling  the  rod  with  more  degrees  of  freedom  it  can  be 


5^ 


shown  that  the  higher  resonant  frequencies  of  the  rod  will 
be  lower  if  the  shaker  and  Impedance  head  are  attached. 

The  magnitude  and  phase  of  the  measured  admittances  are 
affected  appreciably  by  the  mass  of  the  shaker  and  impedance 
head  when  Gd  =  iO*. 


Ct  4*  \ "  I 

iu>xTns) 


(4.21) 

(4.22) 


By  comparing  the  measured  admittances  for  the  total 
ensemble  of  the  rod,  rubber  band,  shaker  and  impedance  head 
to  admittances  calculated  for  a  40-inch  free  rod,  the  effects 
of  the  suspension  system,  shaker  and  Impedance  head  on  the 
measured  admittances  are  evident.  See  Pigures  4.6  through 
4.9. 

The  calculated  admittances  are  generated  by  solving 
equation  4.23  for  a  free  rod. 


E  : 

&  t  C 
3  x2 

St1  u 

(4.23) 

3^Co)  _  ^ 

a%  ~ u- 

(4.24) 

Su.(L')  _  0 
3x  • 

(4.25) 

r  Schematic  of 

- v; 

£  a  Free  Rod 

^  L  * 

Figure  4.5 

L0G-MA6C IN/SEC-LB)  AND  PHASE C RAD  I ANS ] 
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ADMI TTANCEC 1 1 3  MEASURED  FOR  40  INCH  ROD 
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Figure  4.6 

Note:  On  all  admittance  plots  of  the  type  shown  above, 
the  broken  line  is  phase  and  the  solid  line  is  magnitude. 


LOG-MAGC IN/SEC-LB3  AND  PHASECRADI ANS) 
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Figure  4.7 


LOG-MAGC IN/SEC-LB)  AND  PHASECRADI ANS) 
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PAXXINI  « 


Figure  4.8 


LOG-MAGC IN/SEC-LB)  AND  PHASECRADI ANSJ 
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Figure  4.9 


59 


u  Is  the  axial  displacement  of  the  rod  and  Is  a  function  of 
the  location,  x,  and  time,  t. 


£  -  Young's  Modulus 


29.8  x 10 


6 


psl 


A  -  Cross  Sectional  Area 
p  -  Density 
L  -  Length 


.198  Inf 

lb. 


.284 
40.  In. 


in? 


Properties  of  a  40-Inch 
Steel  Rod 
Table  4.1 


The  solutions  to  equation  4.23  Incorporating  the  boundary 
conditions  In  equations  4.24  and  4.25  are  given  In  equation 
4.26. 

0"(x)  =  COS  (4.26) 


The  admittances  of  the  free  rod  are  constructed  In  a 
manner  similar  to  the  method  used  In  section  2.1. 


Y(<0)  -  f  Wn  fc)  4>X*l) 

ml  h>  ^hlto)w  +  i(uf-u>J) 
m.  =  ^  pAcos**ffl*dx. 


njrf 

=  .02  u)n 


(4.27) 

(4.28) 

(4.29) 


(4.30) 


6o 


The  admittances  of  the  free  rod  are  calculated  for  242 
frequencies  from  50  c.p.s.  to  5000  c.p.s.  These  are  the 
same  frequencies  at  which  the  measured  admittances  are  sam¬ 
pled  and  scaled*  The  damping  is  equivalent  to  1%  of  critical 
damping. 


61 


4.3  Predicting  the  Admittances.  Resonant  Frequencies  and 

Mode  Shapes  for  an  80-Inoh  Free  Rod 


The  mutual  and  cross  admittances  for  a  40-inch  free  rod 
are  measured  as  shown  In  Figure  4.10.  The  mutual  admittance 
Is  measured  using  the  accelerometer  in  the  impedance  head. 
The  cross  admittance  is  measured  using  an  accelerometer 
attached  to  the  top  end  of  the  rod. 


T 


ACCELEROMETER 


a.) 


IMPbOA.NCt 

HEAD 

SHAKER 


40-INCh 
^  Rod 


b.) 


A 


A 


T 


&•)  mutual  admittance 
b.')  CROSS  admittance 


F  III  A 

VC1  '  1  vc 

Schematic  of  Measuring  the  Admittances  for  a  40-Inch  Rod 

Figure  4.10 

Having  measured  the  mutual  and  cross  admittances  is 
sufficient  to  construct  two  identical, 2*2,  substructure 
admittance  matrices.  The  mutual  admittance  of  one  end  of  the 
rod  is  the  same  as  the  mutual  admittance  of  the  other  end; 
the  cross  admittances  are  the  same  due  to  symmetry  in  the 
admittance  matrix.  Assuming  that  admittance  measurements 
made  on  one  40-inch  rod  should  not,  in  principle,  differ 
from  admittance  measurements  made  on  a  second  40-inch  rod, 
two  identical  substructure  admittance  matrices  are  con¬ 
structed  from  admittance  measurements  made  on  one  40-inch 


rod 
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The  two  admittance  matrices  for  the  40-inch  rod  sub¬ 
structures  are  inverted  to  substructure  impedance  matrices. 
The  substructure  impedance  matrices  are  assembled,  as  shown 
in  section  2.2,  to  yield  the  impedance  matrix  of  the  compos¬ 
ite  structure,  an  80-inch  free  rod. 


©f__ 


4o-i  ucu 


Schematic  of  Assembling  an  80-Inch  Rod  from  Two  40-Inch  Rods 


Figure  4.11 


The  impedance  matrix  for  the  80-inch  rod  is  inverted  to 
the  composite  structure  admittance  matrix.  The  four  unique 
elements  of  the  composite  structure  admittance  matrix  are 
plotted  as  a  function  of  the  frequency  of  excitation  in 
Figures  4.12  through  4.15.  Figures  4.16  through  4.19  are 
calculated  admittances  for  an  80-inch  free  rod.  These  admit¬ 
tances  are  calculated  using  equation  4.27.  The  damping  in 
the  calculated  admittances  is  1%  of  critical  damping. 

The  admittance  matrix  of  the  composite  structure  is 
used  to  estimate  resonant  frequencies  and  mode  shapes  of  the 
composite  structure  as  shown  in  section  2.4.  The  predicted 
mode  shapes  for  the  80-inch  rod  that  have  corresponding 
resonant  frequencies  that  coincide  with  resonant  frequencies 
of  the  substructures  are  not  accurate.  See  Table  4.2. 


LOG-MAGC IN/SEC-LB)  AND  PHASECRADI ANSD 


Figure  4.12 


LOG-MAGC IN/SEC-LBD  AND  PHASE C RADIANS} 
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PANNlMii  .4 


Figure  4.13 


LOG-MAGC IN/SEC-LB3  AND  PHASE C RAD I  ANSI 
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Figure  4.14 


LOG-MAGC IN/SEC-LBD  AND  PHASECRAOI ANS) 
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Figure  4.15 


log-magc IN/SEC-LB)  and  phasecradi ansd 


67 


Figure  4.16 


LOG-MAGC IN/SEC-LBD  ANO  PHASE C RAO  I ANS ) 
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Figure  4.17 


LOG-MAGC IN/SEC-LB)  AND  PHASECRADI ANSD 
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Figure  4.18 


LOG-MAGC IN/SEC-LBD  AND  PHASEC RAD  I ANS ) 
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PAMS  IS,  4 


Figure  4.19 


Calculated  and  Predicted  Mode  Shapes 
and  Resonant  Frequencies 
Table  4.2 
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Mode 

Shape 

1 


2 


3 


4 


Location 

1 

2 

3 

Resonant 

Admittance 

11 

21 

31 

Frequency 
(c.p.s. ) 

Predicted 

1.0 

* 

.053 

-  .998 

1134. 

Calculated 

1.0 

0.0 

-  1.0 

1257. 

Predicted 

1.0 

-.238 

.142 

2422** 

Calculated 

1.0 

-  1.0 

1.0 

2515. 

Predicted 

.996 

* 

.025 

-  1.0 

3527. 

Calculated 

1.0 

0.0 

-  1.0 

3772. 

Predicted 

1.0 

-  .012 

.012 

4815* 

Calculated 

1.0 

-  1.0 

1.0 

5030. 

*  No  local  maxima  appear  in  this  element  of  the  predicted 
admittance  matrix  at  or  near  this  frequency.  See  section 
2.4. 


**  These  frequencies  are  very  close  to  resonant  frequencies 
of  the  subs true tures. 
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Equations  4.15  through  4.17  show  that  the  measured 
admittances  of  the  40-inch  rod  substructure  are  affected  by 
the  damping  of  the  suspension  system  at  the  first  resonant 
frequency  of  the  free  rod.  In  addition  equation  4.22  shows 

that  the  mass  of  the  shaker  and  impedance  head  also  distorts 

£ 

the  admittances  at  the  frequency  .  Consequently,  the 
admittances  predicted  from  admittance  measurements  made  near 
the  resonant  frequencies  of  the  substructures  are  not  accur¬ 
ate.  Therefore,  if  a  resonant  frequency  of  a  substructure 
is  close  to  a  resonant  frequency  of  the  composite  structure, 
one  should  not  expect  to  predict  accurately  the  corresponding 
mode  shape  for  the  composite  structure. 

In  Table  4.2  the  calculated  mode  shapes,  $(X)  ,  are  de¬ 
fined  in  equation  4.26  with  L  equal  to  80  inches.  The 
calculated  resonant  frequencies  are  defined  as  follows s 


(4.31) 


The  predicted  mode  shapes  and  resonant  frequencies  are  calcu¬ 
lated  as  shown  in  section  2.4. 
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4.4  The  Effects  of  the  Suspension  System  on  the  Measured 

Admittances  of  a  32-Inch  free  Bean 

The  Beam  is  suspended  on  three  lengths  of  £-inoh  shook 
cord  to  simulate  a  free  beam.  The  shock  oord  hangs  parallel 
to  the  T  axis  in  figure  4,33*  See  also  figures  4.20  and 
4.21.  The  beam  is  suspended  on  the  shock  oord  such  that  the 
majority  of  the  weight  of  the  beam  is  balanced  by  the  reac¬ 
tion  force  in  the  length  of  oord  attached  to  the  beam  nearest 
to  the  point  of  forced  excitation.  Consequently,  the  suspen¬ 
sion  is  changed  every  time  the  shaker  and  impedance  head  are 
attached  at  a  new  location. 

figure  4.23  is  an  analog  plot  of  the  mutual  admittance 
of  looatlon  2  on  the  32-lnch  beam  using  a  suspension  system 
similar  to  the  one  shown  in  figure  4.20.  figure  4.24  is  an 
analog  plot  of  the  same  admittance  measured  using  a  suspen¬ 
sion  system  that  supports  the  majority  of  the  weight  of  the 
beam  at  locations  not  near  the  point  of  forced  excitation. 

Two  lightweight  bolts,  1.5  Inches  long,  are  put  into 
the  beam  some  distance  away  from  and  on  both  sides  of  the 
shaker  and  Impedance  head  attachment  point.  One  length  of 
shock  cord  is  attached  to  the  upper  flange  of  the  beam  just 
above  the  Impedance  head  and  shaker,  and  the  other  two 
lengths  of  cord  are  attached  to  the  ends  of  the  bolts, 
figure  4.25  shows  an  X,Y  plane  projection  of  the  suspension 
system  with  the  shaker  and  lmpedanoe  head  attached  to  the 


beam 
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Suspension  System  Used  to  Measure  the 
Admittances  of  a  32-Inch  Free  Beam 


Figure  4.20 
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Measuring  a  Mutual  Admittance 
of  a  32-Inch  Free  Beam 


76 


77 


$199030 


DNVmW!  W  VNVN 


f  ltQUENCY  { Hi) 


78 


S139D30 


3:>NV  13dWI  1 V3INVH  ")}w 


FREQUENCY  (Hi) 


79 


Schematic  of  32-Inch  Beam  Suspension  System 
Figure  4.25 

R*L  is  the  deadweight  reaction  force  in  the  shock  cord 
attached  to  the  beam  near  the  point  of  forced  excitation. 

Rg  is  the  sum  of  the  deadweight  reaction  forces  in  the  shock 
cords  attached  to  the  two  bolts.  and  Wg  are  the  weights 
of  the  shaker  with  impedance  head  and  the  beam,  respectively. 
The  location  of  the  center  of  gravity  for  the  beam,  shaker 
and  impedance  head  is  Xt. 

Using  force  equilibrium  in  the  Y  direction,  equation 

4.32,  and  moment  equilibrium  in  the  Z  direction,  equation 

4.33,  the  deadweight  reaction  forces,  R^  and  R2,  are  deter¬ 
mined. 


(R2  +  k,-v;-v2)  =0. 


.  ( .&■+  Xi) 

CK+V*)  (V Jit) 


(4.34) 
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/  (A  ~  jO 
R2=  (vi+Wz  )(A,+jez) 


(4.35) 


Because  E2  is  the  sum  of  the  two  deadweight  reaotion 
forces  in  the  shock  cords  attached  to  the  bolts  and  each  of 
the  two  deadweight  reaction  forces  at  the  bolts  must  be  pos¬ 
itive  for  the  suspension  system  to  be  statically  stable,  R2 
is  the  largest  possible  force  in  either  of  the  two  shock 
cords  attached  to  the  bolts.  In  general,  the  deadweight 
reaction  forces  in  the  cords  attached  to  the  bolts  are  dif¬ 
ferent. 

During  the  forced  excitation  of  the  beam  the  force  in 
the  shock  cord  changes  very  little,  assuming  that  the  motion 
of  the  shock  cord  is  small.  Assuming  the  angle  between  the 
shock  cord  and  the  vertioal,  0  ,  is  a  function  of  the  dis¬ 
placement  of  the  beam  where  the  cord  is  attached,  one  can 
solve  for  the  force  exerted  in  the  X  direction  by  the  cord 

on  the  beam,  F_. 

it 


Shock: 

Cord 


Schematic  of  the  Reaction  Forces  in  the  Suspension  System 

Figure  4.26 
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(4.36) 


Fa  is  proportional  to  the  deadweight  reaotlon  foroe,  R,  and 
the  looal  displacement  of  the  beam,  €.  • 

To  simulate  a  free  beam  it  is  best  to  make  the  dead¬ 
weight  reaction  forces  as  small  as  possible.  The  suspension 
system  shown  in  Figure  4.21  is  an  attempt  to  minimize  Eg  and 
make  R^^W^+W^.  Therefore,  the  reaction  forces,  F^,  are 
small  at  locations  away  from  the  point  of  excitation  and  the 
reaotlon  force  at  the  point  of  excitation  is,  hopefully, 
much  smaller  than  the  force  applied  by  the  shaker. 

If  Fg  is  the  foroe  applied  to  the  beam  in  the  X  direc¬ 
tion  by  the  shaker  and  F^g  is  the  reaction  force  in  the  X 
direction  exerted  on  the  beam  by  the  shock  cord  attached  at 
the  point  of  excitation,  one  can  solve  for  the  measured 
mutual  admittance. 


(4.37) 


F  =  - 

rRs  x 


(4.38) 


£5  is  the  local  displacement  of  the  beam  at  the  point  of 

excitation  and  Rg  is  the  deadweight  reaction  force  in  the 

shock  cord  attached  to  the  beam  at  the  point  of  excitation. 

If  R  =  W  +  W.,  a  bound  for  the  effect  of  the  shock  cord 
S  12* 

on  the  measured  mutual  admittance  is  established. 


(4.39) 


Wl£Wj)’L  |Ym+' 

AUH  J  i  +YUC 


For  higher  frequencies  of  excitation  the  effect  of  the  dead¬ 
weight  reaction  force,  Rg,  is  very  small  except  at  resonant 
frequencies  of  the  beam. where  1®  also  very  small. 

“true 

Consequently,  the  suspension  system  affects  the  admittance 
measurements  more  at  or  near  the  lower  resonant  frequencies 
of  the  32-inch  beam? 

The  frequencies  of  local  maxima  in  the  measured  admit¬ 
tances  that  have  an  accompanying  phase  shift  of  7T.  radians 
are  compared  to  calculated  resonant  frequencies  for  the  32- 
lnch  free  beam. 


Resonant 

Frequency 

1 

2 


0 

Calculated 
224.  c.p.s. 
628.  c.p.s. 


Measured 
206.  c.p.s. 
549.  c.p.s. 


Comparison  of  Measured  and  Calculated  Resonant 
Frequencies  for  a  32-Inch  Free  Beam 
Table  4.3 


*  See  reference  8  for  the  calculated  admittances  of  a  free 


beam 
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Figure  4.27 


LOG-MAGC IN/SEC-LBD  AND  PHASECRADI ANS) 
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F  ANN  I  MO  004 


Figure  4.28 


LOG-MAGC IN/SEC-LB)  AND  PHASECRADI ANS) 
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r*NMtNQ  004 


Figure  4.29 


LOG-MAGC IN/SEC-LBD  AND  PHASECRADI ANS3 
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fAHNinO  on« 


Figure  4.30 


LOG-MAGC IN/SEC-LB]  AND  PHASECRADI ANSD 
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Figure  4.31 
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4.5  Predicting  the  Admittances,  Resonant  Frequencies  and 

Mode  Shapes  for  a  42-Inch  Oomposlf  Beam 


The  admittances  of  a  42-inch  composite  beam  are  pre- 
dloted  from  admittance  measurements  made  on  a  32-inoh  beam. 
The  32-lnch  beam  has  a  channel  cross-section  and  Is  made  of 
steel. 


Half  Soale  Gross  Section  of  the  Channel 
Figure  4.32 


Iyy  -  Moment  of  Area 
A  -  Oross-Sectional  Area 
E  -  Young's  Modulus 
P  -  Density 
L  -  Length 


.079  in* 
.751  in? 
29.8  *  10* 

*284  l^f 

32.0  in. 


psi 


Physical  Properties  of  a  32-Inch  Beam 

Table  4.4 


Admittance  measurements  are  made  at  three  locations  on 
the  32-lnch  beam  in  the  X  direction  as  shown  in  Figure  4.33 
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In  general,  nine  admittances,  three  mutual  admittances  and 
six  cross  admittances,  are  needed  to  construct  a  3x3  sub¬ 
structure  admittance  matrix.  Because  the  admittance  matrix 
is  symmetric,  one  needs  to  measure  only  six  admittances. 


The  Locations  of  the  Measured  Admittances  for  a  32-Inch 

Beam 

Figure  4.33 

The  mutual  admittance  at  location  1  is  the  same  as  the 
mutual  admittance  at  location  3,  so  the  number  of  admittances 
that  need  to  be  measured  is  reduced  to  five.  The  five 
measured  and  sampled  admittances  are  in  Figures  4.27  through 
4.31.  Each  admittance  is  sampled  for  249  frequencies  be¬ 
tween  50  c.p.s.  and  1000  c.p.s. 

Two  identical  substructure  admittance  matrices  are  con¬ 
structed  from  the  measured  admittances  of  the  32-inch  beam. 
The  substructure  admittance  matrices  are  inverted  to  sub¬ 
structure  impedance  matrices  which  are  assembled  into  a 
composite  structure  impedance  matrix  as  shown  in  section  2.2. 
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The  composite  structure  is  shown  schematically  in  Figure 


Assembling  a  42-Inch  Composite  Beam 
from  Two  32-Inch  Free  Beams 
Figure  4.34 

The  composite  structure  is  made  of  two  beams  pinned 
together  at  two  points.  At  these  two  points  of  attachment 
the  beams  have  compatible  X  displacements;  however,  the  beams 
rotate  independently  about  the  T  axis  at  these  two  points. 
Compatibility  of  displacement  in  the  Z  direction  at  the 
attachment  points  is  not  necessary,  since  the  composite 
structure  has  no  degrees  of  freedom  in  the  Z  direction. 

The  composite  structure  impedance  matrix  is  inverted  to 
a  composite  structure  admittance  matrix  and  the  mode  shapes 
and  resonant  frequencies  of  the  composite  structure  are 
estimated  as  shown  in  section  2.4.  The  unique  elements  of 
the  composite  structure  admittance  matrix  appear  in  Figures 
4.35  through  4.40. 

The  predicted  resonant  frequencies  and  mode  shapes  for 
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LOG-FREQ  CCPS) 


Figure  4.35 
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Figure  4.36 
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LOG-FREO  (CPS!) 


Figure  4.37 
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Figure  4.38 


LOG-MAGC IN/SEC-LBD  AND  PHASECRADI ANSD 
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Figure  4.39 


LOG-MAGC IN/SEC-LBD  AND  PHASECRADI ANS) 
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Figure  4.40 
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the  composite  structure  are  compared  to  resonant  frequencies 
and  mode  shapes  calculated  for  the  composite  structure. 

The  calculated  data  Is  generated  by  a  finite-element  com¬ 
puter  program,  STRUDL  II,  developed  by  the  MIT  Oivil  Engi¬ 
neering  Department^  ^  A  listing  of  the  predicted  and  cal¬ 
culated  data  appears  In  Table  4.5. 

Two  of  the  predicted  resonant  frequencies  of  the  com¬ 
posite  structure  are  very  olose  to  resonant  frequencies  of 
the  substructures.  As  shown  In  section  4.3  the  predicted 
mode  shapes  for  these  two  resonant  frequencies  of  the  com¬ 
posite  structure  should  not  be  accurate  due  to  the  distortion 
of  the  measured  admittances  at  the  substructure  resonant 
frequencies.  The  two  resonant  frequencies  of  the  composite 
structure  that  are  near  resonant  frequencies  of  the  sub¬ 
structures  correspond  to  the  second  and  fourth  mode  shapes 
In  Table  4.5. 

Equation  4.39  shows  that  the  effect  of  the  suspension 
system  on  the  admittance  measurements  should  decrease  with 
Increasing  frequency  of  excitation.  If  the  suspension 
system  Is  responsible  for  the  Inaccuracy  of  the  second  and 
fourth  predicted  mode  shapes,  the  second  mode  shape  should 
be  less  accurate  than  the  fourth  mode  shape.  The  data  In 
Table  4.5  shows  that  the  second  predicted  mode  shape  at  202 
c.p.s.  Is  less  accurate  than  the  fourth  predicted  mode  shape 
at  549  c.p.s. 

The  calculated  resonant  frequencies  for  the  composite 
structure  are,  In  general,  higher  than  the  resonant  frequen- 
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Calculated  and  Predicted  Mode  Shapes  and  Resonant 
Frequencies  for  a  42-Inch  Composite  Beam 

Table  4.5 


Location 

1 

2 

3 

4 

Resonant 

Admittance* 

11 

21 

31 

41 

Frequency 
(c.p.s. ) 

Mode 

Shape 

1  Predicted 

1.0 

-  .052 

-.052 

1.0 

130. 

Calculated 

1.0 

-.059 

-.059 

1.0 

147. 

2  Predicted 

-1.0 

-.828 

.669 

-  .715 

*# 

202. 

Calculated 

-1.0 

.676 

-.676 

1.0 

205. 

3  Predicted 

-  .998 

.454 

-.453 

1.0 

305. 

Calculated 

-1.0 

.496 

-.49  6 

1.0 

415. 

4  Predicted 

-.893 

.839 

.578 

-1.0 

#* 

549. 

Calculated 

-.986 

1.0 

1.0 

-.986 

567. 

5  Predicted 

.996 

-.325 

-.330 

1.0 

851. 

Calculated 

1.0 

.178 

.178 

1.0 

1022. 

6  Predicted 

-  .971 

.391 

.406 

-  1.0 

870. 

Calculated 

-1.0 

.294 

-.294 

1.0 

1037. 

*  The  admittances  are 

used  only 

for  the 

predicted 

data. 

These  frequencies  are  near  resonant  frequencies  of  the 
32-inch  beam  substructure. 
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cles  predicted  by  the  admittance  data.  The  composite  struc¬ 
ture  Is  modeled  using  the  assumed  displacement  method  which 
should  yield  calculated  resonant  frequencies  that  are  upper 
bounds  for  the  resonant  frequencies  of  the  real  structure. 
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4.6  Conclusions 


The  resonant  frequencies  and  mode  shapes  of  a  composite 
structure  can  be  predicted  using  admittance  measurements 
made  on  the  constituent  substructures  If  certain  conditions 
are  met.  First,  the  shaker  must  be  able  to  produce  enough 
force  to  measure  the  admittances  of  the  substructures  over 
the  frequency  range  of  interest.  Second,  the  composite  must 
not  have  resonant  frequencies  that  are  very  close  together. 
Third,  the  suspension  system  must  not  alter  the  admittance 
measurements  on  the  substructures.  Finally,  the  composite 
structure  must  not  have  resonant  frequencies  that  coincide 
with  resonant  frequencies  of  the  constituent  substructures. 
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Appendix  A 


The  method  of  matrix  inversion  used  in  the  computer 
program  ADMIT  is  "triangular  decomposition. This  method 
of  matrix  inversion  entails  decomposing  a  square  symmetric 
matrix  into  the  matrix  product  of  a  matrix  [t>l  and  its 

transpose . 


(y]  =  [D][Df  (a-1) 

The  matrix  [D]  is  constructed  such  that  all  elements  to  the 
right  of  the  principal  diagonal  are  zero.  There  are  the  same 
number  of  unique  elements  in  to]  as  there  are  unique  ele¬ 
ments  in  M  .  The  elements  of  Co]  are  solved  for  using 
equations  A. 2  through  A. 5. 


(A.2) 
(A. 3) 


n-» 


dvntx  “  c)  p]/ VY\>V\ 


(A.4) 

(A. 5) 


The  inverse  of  M  is  written  as  a  product  of  the  in¬ 
verse  of  [0]  and  the  transpose  of  the  inverse  of  to]. 


Cy  r'  -  [d-TV 


(A.  6) 


[E] -  [0] 


(A. 7) 
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The  elements  of  the  matrix  It  J  are  solved  for  using  the 
definition  of  an  Inverse  matrix,  equation  A. 8. 

[D][0r=  [D)[E}=  [l]  U.8) 

Since  all  of  the  elements  of  Coi  to  the  right  of  its 
principal  diagonal  are  zero,  it  can  be  shown  that  all  of  the 
elements  of  (El  to  the  right  of  the  principal  diagonal  are 
zero.  Given  all  the  elements  of  [o}  ,  one  can  solve  for  the 
elements  in  [E]  using  equations  A*9  and  A*10. 


rn>n. 


(A. 9) 
(A. 10) 
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Appendix  B 

Computer  Program  Listings 
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